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Abstract 

The temperature dependence of the I-V characteristics of many single- 
electron tunneling devices enable thermometer operation of these systems. 
We investigate two normal conducting kinds of them, (a) a, single junction 
in a high-impedance environment, and (b) a double junction. The char- 
acteristics of both devices show a crossover from Coulomb blockade at low 
temperatures to ohmic behavior at high temperatures. The related differen- 
tial conductivity dip allows the determination of the junctions temperature. 
Both configurations (a) and (b) are expected to work at least within the range 
0.5 < PEc < 10, where Eq is the Coulomb energy of the system under investi- 
gation. We present an analytical solution for both low- and high-temperature 
case of (a) and (b) as well as numerical results and their fit, including the effect 

of co-tunneling in case of a double junction. 
PACS 73.40G, 73.40R 
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I. INTRODUCTION 



Sensoring turns out to be one of the possible applications of single-electron tunneling 
(SET) devices. After the establishment of the electrometer that is one of the basic 
devices of SET, radiation detectors and resistor comparators were studied. Nowadays, 
SET-devices as thermometer are brought into the scope of investigation [^. 

For this operation the known temperature dependence of the current through ultrasmall 
tunnel junctions is exploited. As the single junction without any environment, i.e. the 
voltage biased tunnel junction, does not show a reasonable influence of the temperature 
on its characteristic, the basic thermometer devices are (a) the normal conducting single 
junction in a high-impedance environment and (b) the normal conducting double junction. 
Studies on both ('aj |^ and fbj are done. This paper is restricted to (a) the infinite high- 
impedant ohmic environment (-Re ^ oo) and (b) the symmetric double junction {Ri = R2, 
Ci = C2). On the other hand, it extends because of spotting the thermometer operation 
as well as due to a more rigid theoretical treatment. 

The suppression of thermal fluctuations in SET junctions is expressed by the condition 
Eq ^ k-gT, where Eq is the Coulomb energy of the considered system. If this condition 
holds both (a) and (b) show a Coulomb blockade 0- This is the low-temperature case 
of our consideration. At high temperatures, Eq <^ k^T, Coulomb effects are covered by 
thermal fluctuations and an ohmic current-voltage characteristic remains. 

The thermal dependence is rather discussed in terms of the differential conductivity G{V) 
of the junction(s) than the current. According to the Coulomb blockade G{V) shows a dip 
at \^ = (Fig. [l|, ^. This dip is well-pronounced at low temperatures and is smeared out 
with increasing temperatures. Therefore, the shape of the dip allows the determination of 
the temperature of the junctions. This is the thermometer operation. 

There are different measures of the mentioned dip that may be used for evaluation. In 
this paper the depth of the dip, i.e. G(y = 0), is used only. It is preferred to the width of 
the dip because a precise measurement of the latter at higher temperatures would require 
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improper high voltages across the junction(s). Additionally, the depth of the dip is easier 
accessible in theory. 



II. SINGLE JUNCTION 



A single junction possesses the Coulomb energy Ec = e^/(2C), where C is the junction 
capacitance. The current via this junction in a high-impedance ohmic environment {Re 
oo) is given by 
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in terms of the junction resistance R and the temperature (3 = l/{kBT). To calculate this 
integral the first term of it is represented as 
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This representation breaks down at /? = 0. Furthermore, the convergence of the sum is poor 
in the vicinity of this value. Therefore, this domain is excluded now but studied in detail 
later. Due to the sum, (j^) simplifies to two integrals of the type 
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with erfc(a;) = 1 — erf(x). The corresponding calculation yields 



(3) 



Ec ^ r„/eV^ 



'(V) - ^ E K 







Ec, ' 



n 



^ Ec' 



(4) 



with 



f(x, n) 



-(3Ec{x-lf/i 



V^El 



(5) 



-(2n - X + 1) exp[pEc{2n -x + l)V4]erfc((2n - x + I) 
■{2n + X+1) exp[/?Ec(2n + x+ l)V4]erfc((2 n + X + 1) 



\(3Ec 



ipEc 



Obviously, the current is the difference of two contributions depending on each other by the 
substitution V ^ —V. The determination of the minimal differential conductivity yields in 
terms of 
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If l3Ec > 1 holds the asymptotic expansion 
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can be used to transform into 



it V vr (P-C/c) 
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using the Bernoulli numbers B2m- The term corresponding to m = 1 yields the behavior of 
(p!0|) for /5i?c ^ 1- This is the short-dashed curve in Fig. 0. 
In case of 13 Ec <^ 1 the evaluation of (0) uses 

e 1 



/3g/2 



(11) 



1-e-^^ /3 sinh/3£/2' 

Therefore, the integral in (|I]) includes besides the Gaussian an other peak of the shape 
x/ sinhx. Because the width of the Gaussian and the x/ sinhx-peak are proportional to 
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{/3Ec)~^^'^ and {PEc)~^, respectively. Hence, the latter peak becomes broad and fiat com- 
pared to the Gaussian in the limit P ^ 0. Therefore, it turns out to be sufficient to use the 
approximation 

3e/2 , , 

sinh /3£/2 ^ ' 

in this limit and an integral of the Gaussian remains of (|l|). It results in 

liy) = ^(l - e-'^'^^) exp [ - ^(Ec - 2e V)\ (13) 
for the current via the junction and gives 



(14) 



for the minimal differential conductivity of the junction. This formula is represented in 
Fig. by means of a long-dashed curve. 

According to Fig. ^, an evaluation of (|I]) at (HEq ^ 1 is desired, but hardly to achieve 
analytically. Therefore, a simple fit to the numerical data was performed, resulting in 

G(y = 0) = ^ (l - - arctan(0.5/3Ec)) . (15) 

This formula corresponds to the dotted line in Fig. The fit may be enhanced by introducing 
a polynomial into the argument of arctan. In the given simple form (|15D enables an inversion 
to 
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/5Ec = — tan[-(^l-G(\/ = 0)i?JJ, (16) 
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a simple thermometer function. 



III. DOUBLE JUNCTION 

In case of the symmetric double junction the Coulomb energy is given by Eq = e^/ (4C), 
where C is the junction capacitance of each junction. Due to several states of the central 
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electrode of the double junction, that are determined by the number n of extra charges on 
it, the current I{V) turns out to be a sum over contributions of each state : 



^(^) = {^li.^) - h{n) (T{n) = [r2{n) - kin) 



E 



(17) 



In general, both the rates ri 2{n), h,2{n) for both directions of the first and second junction 
and the occupation probability a{n) are voltage dependent. The voltage dependence of the 
rates is given in terms of the energy gain per tunnel event [|T^ 
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where R is the junction resistance. 
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According to the Likharev-Averin equation the stationary occupation probabilities 
(T{n) are governed by 

a{n + 1) ri(n) + l2{n) 



a{n) r2(n + l) + /i(n + l)' 
In case of n ^ 1 this simplifies to 
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Comparing with a Gaussian normal probability function with mean value and width s 
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the occupation probability can be approximated by this Gaussian, if 

1 



(23) 



is chosen. This approximation is in the limit /3 — exact. However, according to Fig. |^ 
it seems to be reasonable, if (5Ec < 10 holds. The main result of this consideration is the 



weak voltage dependence of the occupation probabilities a{n). Neglecting this dependence, 
the differential conductivity of the studied double junction might be expressed as 

Giy) = ^j^ = T.9{nMn), (24) 

where 

9{n) = (25) 

is used with regard to (|T7|) . Investigating the minimal differential conductivity of the double 
junction, i.e. V = 0, the calculation yields 

. _ 1 sinh/3Ec + /3Eccosh/3Ec 

g[n)\y^, - 2i? r cosh f3Ec - cosh 2n(3Ec ^ ^ 

[2n(3Ec sinh 2n(3Ec - f]Ec sinh (3Ec] sinh (3Ec 1 
[cosh (3Ec - cosh 2n(3Ecf J ' 

In case of n = this result simplifies to 

11-effc+^^e^ 
^(°)lv^=o = R (e^i^c - 1)^ • (27) 

Further simplification is achieved in the case of f3Ec ^ 1, and f3Ec ^ 1 where 

(3Ec « 1 : ^(0) - ^ (l - ^) (28) 



/?Ec » 1:^(0)^^ e-^^c 



holds, respectively. In both cases further evaluation is possible owing to (121) and the nor- 
malization of cr{n), 

E^H = 1- (29) 

n 

As g{n) turns out to be independent of n up to second order in PEq the approximation 
g{n) ^ g{0) can be exploited in the first case {PEq <^ 1). This results in 

G{V = 0) ^ g{0) Y: <n) = g{0) = ^ (l - ^) • (30) 
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This formula is represented by the long-dashed curve in Fig. ^. In the opposite limit, i.e. 
PEc ^ 1, o'(n) is sharply peaked at n = (see Fig. ^(c)) and the evaluation yields 

GiV = 0) ^ ^(O)a(O) = g{0) = ^e"^^-. (31) 

This result is shown by the short-dashed curve in Fig. ^. 

Again, the missing analytical approximation at (3Ec is substituted by the following fit : 

^(^ = 0) = - tanh(0.3/3Ec)] (32) 
ZK 

f3Ec = — artanh[l - 2G{V = 0)R] 
0.3 

This fit is given by the dotted curve in Fig. ^. The second line of ( |32[) enables the determi- 
nation of the temperature out of the minimal differential conductivity G{V = 0). 



IV. CO TUNNELING 



According to |]TT| the correlated transfer of electrons across both junctions of a double 
junction results in current contributions due to two different processes known as coherent 
and incoherent co-tunneling. The influence of the coherent co-tunneling can be neglected 
in this consideration because it depends on the size of the central electrode which is a design 
parameter. Hence, it should not be a problem to suppress this influence by means of a larger 
island between the junctions. 

On the other hand, the current contribution of incoherent co-tunneling is independent 



of this size. It is given by |[TT 



The evaluation in case of the symmetric double junction (i?i = R2 = R, AEi\v=o = 
AE2\v=o = Eq) at = yields for the contribution to the differential conductance 
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using the abbreviation Rq = h/e^. In order to calculate the total differential conductance at 
V = 0, this term has to be added to the results of Sec. fT|. Hence, the thermometric use of the 
double junction at higher temperatures (/3-Ec 1) requires high tunnel resistances. Fig. |^ 
shows the high-temperature approximation of the conductance G{V = 0) for different tunnel 
resistances. As seen from Eq. (P^) co-tunneling becomes negligible for low temperatures 



{f3Ec ^1). In case of a realistic value R ^ IORq incoherent co-tunneling restricts the 
application to PEc > 0.5. If co-tunneling reaches the order of the one-particle process, the 
calculation which bases on a perturbation theory becomes insecure. 

V. DISCUSSION 



There are obvious similarities in the behavior of the devices of Sec. || and Sec. |T| to 



build a simple thermometer on base of SET junctions as far as co-tunneling is excluded. 
In this case the characteristics G{(3Eq) are comparable, but due to co-tunneling that rises 
strongly with growing temperature the working field of the double junction is restricted to 
0.5 < (3Eq < 10 approximately, whereas the single junction is expected to work reasonable 
within 0.1 < I3Eq < 10. Therefore, the single junction is preferred from a theoretician's 
point of view. 

In experimental realization, however, it is much harder to prepare a current biased single 
junction than a voltage biased double junction. Hence, an experimenter would prefer the 
second kind of SET-thermometer. According to these points the cons and pros of both 
configurations are balanced. It could be expected that multi junction arrays are better 
suited at high temperatures {PEq < 1) because of suppressed co-tunneling. 

Practical exploitation of SET thermometers requires the knowledge of the junction pa- 
rameters, i.e. the junction capacitance and resistance. These parameters might be deter- 
mined by measuring the I{V) characteristic of the system at low temperatures jSEc ^ 1. 
After calibration the thermometer is ready to be operated at /SEq ~ 1. Therefore, the 
thermometers might be used in the mK-region and their field of application is rather nar- 
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row. Owing to quantum fluctuations the effect will be disturbed for very low temperatures if 
SET-conditions hold (i? > /i/e^). In case of a single junction a finite environment resistance 
will cause changes. 

One possible kind of usage is as on-chip thermometer of other SET experiments. The 
advantage of it is the accurate determination of the junction temperature on the chip. Due 
to (3Ec ~ 1 the technical restrictions of the thermometer junctions are less rigid than those 
of the other SET junctions on chip, where PEc 3> 1 must be fulfilled. The disadvantage 
is the above mentioned calibration that requires even lower temperatures than those of the 
following experiment. 

The presented theoretical calculations allow more detailed treatment of the thermal 
behavior of SET junctions. The search for an analytical approximation at (3Ec ^ 1 might 
be successful. An interesting point of the thermal behavior of small tunnel junctions that 
may become important within a more detailed treatment is the electron heating by the 
current [|T^. However, in the standard framework applied here this effect is neglected. 
Eventually it must be taken into account for lower temperatures {f3Ec ^1). In this case 
the theoretical approach basing on thermodynamic equilibrium is questionable. It should 
be noted that in this case the electron temperature exceeds the substrate temperature [IT2 



In Sec. the study of a non-ohmic or non-infinite environment might make sense. The 
investigation of an asymmetry might be of interest in case of the double junction. Therefore, 
this paper gives a first theoretical approach only. 

For useful hints we are indebted to L. S. Kuzmin and S. Vyshenski. This work was 
supported by the Deutsche Forschungsgemeinschaft. 
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FIGURES 

FIG. 1. Differential conductivity of a single junction within a high-impedance environment for 

varying temperature showing the change of the differential conductivity dip at V = 0. The depth 
of the dip is estimated. The curves belong to (3Ec = 0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1.0, 2.0, 5.0, 10, 
20, 50, and 100. 

FIG. 2. Minimal differential conductivity G{V = 0) of a single junction in a high-impedance 
environment. The solid line links the values taken from Fig. 1, the dotted line corresponds to a fit 
formula given in the text. The long and short-dashed curves show analytical approximations for 
pEc < 1 and pEc » 1. 

FIG. 3. Differential conductivity of a symmetric double junction in dependence on the applied 
voltage for temperatures pEc = 0.1, 0.2, 0.5, 1.0, 2.0, 5.0, 10, 20, 50, 100. The peaks at low 
temperatures arise from the Coulomb staircase. 

FIG. 4. Minimal differential conductivity G{V = 0) of a symmetric double junction. Again, 
the solid line belongs to the data of Fig. 3 and the dotted line to a fit. The long and short-dashed 
curves show analytical approximations for ^Eq ^ 1 and PEq S> 1. 

FIG. 5. Actual occupation probability a{n) of n charges on the central electrode of a symmetric 
double junction and a corresponding gaussian probability distribution for varying temperatures (a) 
PEc = 0.1, (6) 1.0, (c) 10. For /3 — the Gaussian becomes exact. 

FIG. 6. Influence of the incoherent co-tunneling on the differential conductance G{V = 0) for 
different values of R/Rq. As in Fig. 4, the dotted line corresponds to the Fit to guide the eye. The 
other curves show the total differential conductance in high-temperature approximation at V = 
for the ratios R/Rq = 1, 10, 100, 1000 from up to down. 
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